Weak localisation in bilayer graphene 
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We have performed the first experimental investigation of quantum interference corrections to the 
conductivity of a bilayer graphene structure. A negative magnetoresistance - a signature of weak 
localisation - is observed at different carrier densities, including the electro-neutrality region. It is 
very different, however, from the weak localisation in conventional two-dimensional systems. We 
show that it is controlled not only by the dephasing time, but also by different elastic processes that 
break the effective time-reversal symmetry and provide intervalley scattering. 
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The discovery that a single layer of carbon atoms 
(graphene [l[) can be manufactured in such a way that 
the carrier density is controlled by a gate voltage [HQ be- 
came a breakthrough in creating a new, high-quality two- 
dimensional (2D) system. It has remarkable differences 
from conventional 2D systems: a linear energy spectrum 
and chirality of charge carriers. These initiated much in- 
terest (mainly theoretical) in exploring graphene's prop- 
erties. One of them is the manifestation of well known 
[H quantum interference corrections to the conductivity. 

Interference of electron waves that are scattered by dis- 
order and form closed trajectories usually decreases the 
conductivity, resulting in the so-called weak localisation 
(WL) of electrons. Perpendicular magnetic field intro- 
duces a phase difference for interfering carriers and de- 
stroys the WL, making measurements of the (positive) 
magnetoconductivity a sensitive tool of quantum inter- 
ference. Carriers in graphene are chiral, that is, they 
have an additional quantum number, pseudospin, arising 
from the fact that their wavefunctions are composed of 
the contributions from two sublattices. If pseudospin is 
conserved, the backscattering of carriers is forbidden and 
weak antilocalisation rather than WL will occur 0, [H, [|| , 
with a positive correction to conductance and a negative 
magnetoconductance. Furthermore, the quantum inter- 
ference in graphene will be sensitive not only to inelas- 
tic scattering of carriers, which breaks the phase of the 
wavefunction, but also to a number of elastic scattering 
mechanisms, especially atomically sharp and topological 
defects which affect the carrier chirality ]B|, 0| . The energy 
spectrum "warping" of the carriers is also of major impor- 
tance for WL in graphene 0. The first experiments on 
WL in single-layer graphene @, H, have indeed shown 
its high sensitivity to the details of scattering processes, 
even full suppression of the interference effect [7|] . 

In this work we study WL in a bilayer graphene struc- 
ture. It has already been understood that adding a 
second layer of graphene changes the system dramati- 
cally: the linear dispersion relation becomes parabolic 
[HI]. In spite of this, bilayer graphene is expected 
to be very different from conventional 2D systems with 
parabolic spectrum, as well as from single-layer graphene 



L ll , 12 . IT3L Ihj ] (as measurements of the quantum Hall 
effect |15l Il6| have already shown). The carriers remain 
chiral, but the Berry phase of the wavefunction (a mea- 
sure of chirality) in a bilayer is 27T and not tt as in a sin- 
gle layer (which means that the pseudospin turns twice 
as quickly in the plane than the momentum, while in a 
single-layer it is aligned with the momentum). As a re- 
sult, there will be no suppression of backscattering, and 
the quantum correction in a bilayer will have the sign of 
conventional WL 

H El- 

Its magnitude, however, will 
still be very sensitive to different elastic processes: WL 
in each of the two graphene valleys can be totally sup- 
pressed by topological defects [E| and warping of the en- 
ergy spectrum (the latter is ex pec ted to be stronger in 

At the same time, 



a bilayer than in a single layer 
intervalley scattering can partially restore WL. 

Our measurements of the positive magnetoconduc- 
tance of a gate-voltage controlled bilayer structure have 
shown the existence of WL at all studied carrier densi- 
ties. The detailed analysis of the results is performed 
using the theory [ijj], which has allowed us to determine 
three physical quantities responsible for the manifesta- 
tion of WL: the phase breaking time and two elastic 
times - the intervalley scattering time Tj and an effective 
time-reversal symmetry breaking time t w . 

The sample of rectangular shape (1.5 /im xl.8 /xm) 
is fabricated by the method of mechanical exfoliation of 
(HOPG) graphite developed in The bilayer flake is 
deposited on a n + Si substrate separated by a 300 nm 
Si02 layer. Using e-beam lithography, two Au/Cr con- 
tacts were made at opposite edges of the flake to allow 
resistance measurements from 40 mK to 4.2 K. The car- 
rier density up to 1.5 X 10 12 cm -2 is varied by the gate 
voltage between the flake and n + Si (in accordance with 
the relation dn/dV g = 7.3 x 10 10 cm~ 2 V determined 
by the thickness of Si02). 

Figure 1(a) shows the dependence of the sample re- 
sistance on gate voltage, with a characteristic peak in 
the region of electro-neutrality, where the type of carrier 
changes from electrons (at positive V g ) to holes (at nega- 
tive V g ). (A small offset in the position of the maximum 
of ~ —3 V, caused by unintentional doping of the flake 
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FIG. 1: (color online) (a) Dependence of the resistance of the 
bilayer flake on gate voltage at 40 mK. (b) Conductance as a 
function of V g at B = 14.3 T. (c) Conductivity fluctuations as 
a function of magnetic field for different temperatures (region 
II). The rectangle defines the range of fields used for WL 
measurements. 



is compensated in Fig. 1 (a).) The carrier mobil- 
ity outside the electro-neutrality region, \V g \ > 5 V, is ~ 
7000 cm 2 V -1 s _1 . The gate voltage dependence of the re- 
sistance in strong magnetic field has shown quantum Hall 
plateaus that indicate that the sample is indeed bilayer 
graphene: the conductance quantisation a xy — ANe 2 /h, 
where N is an integer, occurs at Landau-level filling fac- 
tors v = 1/2,3/2,... [l6j]. An example of such depen- 
dence with a clear v — 1/2-plateau is given in Fig. 1(b). 

In our study of WL we concentrate on a small mag- 
netic field range indicated by the box in Fig. 1(c). Be- 
cause of the mesoscopic nature of the small sample, its 
conductance shows reproducible fluctuations of the order 
of e 2 /h, both as a function of V g and B, decaying with 
increasing temperature, Fig. l(a,c). The presence of the 
fluctuations made the magnetoconductance traces sensi- 
tive to small shifts in V g . Therefore, we have averaged 
the effect of magnetic field over the range AV g = 2 V, 
which is not too large to change significantly the average 
resistance of the sample, but large enough to contain sev- 
eral R (V g ) -fluctuations corresponding to several "sample 
realisations". Five studied V^-regions are shown in Fig. 
1(a). 

Figure 2 illustrates our procedure of averaging. At 
a given temperature, we find the conductivity cr(V g ) at 
several values of B- field, which is changed up to ± 17 
mT with a step of ~ 1 mT. Then the average difference 
Aa(B) = {a(B,V g ) - <r(0, V g )) A v 9 is determined. An 
example of the resulting magnetoconductivity (MC) at 
several temperatures for regions I and II is presented in 
Fig. 3. To decrease further the effect of conductance fluc- 
tuations, the data in region II are averaged with those in 
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FIG. 2: Evolution of conductance fluctuations with small 
magnetic fields applied, for regions I and II. 



region V which is approximately symmetric with respect 
to the electro-neutrality point. (To find a from the mea- 
sured resistance R(V g ), the resistance of the contacts, 
~ 175 Ohm, is taken into account. This value is found 
from the deviation of the plateau in Fig. 1(b) from the 
expected quantised value of 4e 2 /h.) 

To analyse the results we use the following expression 
for the MC due to WL in bilayer graphene 141 ] : 
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Here ip(x) is the digamma function, is the phase 
breaking time, t% is the intervalley scattering time, and 
(r*) _1 = (r.() _1 + (ru,)" 1 , where t w is the intra- valley 
"warping" time combined with the time of chirality 
breaking Equation (1) contains three terms, and its 
main difference from the single-layer result @ is in the 
positive sign of the third term. It is seen that in the case 
of Ti,T* — * oo, i.e. in the absence of the corresponding 
elastic scattering, the MC is reduced to the third term 
and becomes the conventional expression for WL in a 
2D system with two valleys, with a positive MC. On the 
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FIG. 3: (color online) Averaged magnetoconductivity for re- 
gions I and II. Dashed curves are the fits using only the first 
term in Eq. (1), and solid lines are the fits with the first two 
terms. 



contrary, strong warping scattering (t w — > 0) without in- 
tervalley scattering (ji — * oc) suppresses the third term 
while the first two cancel each other - the result is zero 
MC (complete suppression of WL). However, introduc- 
ing some intervalley scattering (t; ~ t^) will make the 
second term smaller than the first, so that there will be 
no full cancellation and the first term will dominate in 
small fields. The result will be a conventional, positive 
MC, although with a suppression factor of 1/2. We will 
now see that it is exactly this situation that describes the 
experimental results. 

Figure 3 shows that the MC in very small fields can 
be indeed satisfactorily approximated by the first term 
in Eq. (1), but with increasing B this approximation 
strongly deviates towards larger values. Inclusion of the 
second (negative) term allows one to get good agreement 
with experiment over the whole field range and obtain 
values of the temperature-dependent dephasing time 
and temperature-independent intervalley time Tj. (Some- 
what better agreement can be achieved by including the 
third term of Eq. (1), from which the value of the warp- 
ing time can be estimated - see below.) 

We start the discussion of the results by noting that 
the conductivity of our sample is always larger than e 2 /h, 
Fig. 1, and therefore the pertubation theory of WL is 
applicable in the whole range of V g . Analysis of the data 




FIG. 4: (color online) (a) Temperature dependence of the 
dephasing rate for regions I (solid line) and II (dashed line) 
fitted with a linear dependence. Inset: Coherence length as 
a function of temperature for the two Vg-regions. (b) Coher- 
ence length at 40 mK (open circles) and intervalley scattering 
length (half-filled circles) for different Vg-regions. (Dotted line 
is a guide to the eye.) c) Dephasing time at 40 mK (open cir- 
cles) and 4.2 K (squares) in comparison with the intervalley 
time (half-filled circles) and momentum relaxation time (filled 
circles), (d) Comparison of the fits with two terms (dotted 
line) and three terms (solid line) in Eq. (f), for the 40 mK 
curve in region I, Fig. 3. 



with the first two terms in Eq. (1) (assuming strong 
warping, t w — > 0) gives the value of the coherence length 
= (Dt^) 1 / 2 and the phase breaking time t^, Fig. 
4(a). The diffusion coefficient is found from the relation 
D = airh 2 /2m*e 2 , where a is the average conductivity 
in each region and to* = 0.033 mn is the effective mass 
of carriers in graphene bilayers 11 1 . 

The manifestation of WL in region I deserves special 
consideration. It is reasonable to assume that when the 
average carrier density is close to zero, the presence of 
disorder splits the system into electron- hole puddles 17 1 
separated by p-n junctions 12, HI- (A rough estimation 
of the carrier density in the puddles, n,p < 3 x 10 cm~ 2 , 
can be obtained from the range \V g \ < 4 V where, ignor- 
ing the fluctuations, the dependence R(V g ) in Fig. 1 
flattens out.) The fact that WL in region I is not sup- 
pressed means that breaking down the system into pud- 
dles does not prevent interference of charge carriers on 
a scale larger than the puddle size. A possible reason 
for this is that the p-n junctions separating the puddles 
are highly transparent (at some angles of incidence) for 



4 



chiral carriers 12|, |l8j, and therefore carrier transitions 
between the puddles can occur in a coherent way. 

Figure 4(a) compares L<p(T) and t^{T) in regions I and 
II. It is seen that the coherence length increases with de- 
creasing T but approaches a finite value at T — ► 0, which 
is also seen as an offset in t~~ 1 (T) . (Our preliminary anal- 
ysis of the <j{B)- fluctuations as universal conductance 
fluctuations has given similar magnitude and behaviour 
of L(f,(T).) There is a natural explanation that the in- 
crease of L$ at low temperatures is limited by the size of 
the sample, L = 1.5 fim. The somewhat smaller satura- 
tion value of L<j, in region I can be treated as an indica- 
tion of the formation of puddles in the electro-neutrality 
region. Charge carriers inside a puddle have to scatter 
many times from its boundary before finding an optimal 
angle to penetrate the p-n junction. This results in an 
effective decrease of their coherence length. 

The increase of the dephasing rate with temperature 
in all regions agrees with the prediction of the theory of 
electron-electron interaction in the diffusive regime [l9| , 
where ^TT p /fi < 1. This regime is indeed realised in 
our sample where the parameter ksTTp/h varies from 
0.001 to 0.1. We have approximated the experimental 
temperature dependence of the dephasing rate at T > 
1 K as r" 1 = phsT In g/(hg), where g = crh/e 2 is the 
dimensionless conductivity. The empirical coefficient (3 
is found to be close to unity: ~1.2 in all regions, Fig. 
4(a). 

Figure 4(b) shows the low-T saturation value of L^. 
It is close to half the length of the sample in all regions 
apart from region I where it gets smaller, which we ten- 
tatively ascribe to the formation of the puddles. Also 
shown is the temperature independent intervalley diffu- 
sion length U = (Dt;) 1 / 2 . This can be interpreted as 
a characteristic distance between defects that are able 
to change significantly the momentum of charge carriers. 
By comparing this distance with the size of the sample, 
one can estimate the number of such defects: ~ 10. 

Figure 4(c) compares, in several V 9 -regions, the phase 
breaking time with the intervalley scattering time n 
and momentum relaxation time t p . One can see that 
70 and 7j are of the same order of magnitude and much 
larger than r p ~ 0.15 ps. An estimation of the warping 
time can be obtained by taking into consideration the 
third term in Eq. (1). It has less effect on the description 
of the MC and does not change much the values of 
and n found earlier. Figure 4 (d) illustrates the degree 
of improvement of the fit with inclusion of the third term. 
By considering its contribution, we obtain the following 
estimation of the warping time in the studied V g -regions: 
t* ~ t w < 0.5 ps, which is close to the estimation [14| of 
the warping time in bilayer graphene at densities ~10 12 
cm~ 2 : t w ~ T p . The small value of the obtained t w is 
consistent with the smallness of the third term and is a 
clear signature of the importance of the warping in the 
studied system. 



It is interesting to compare our values of r, with that 
recently obtained from the analysis of WL in single-layer 
graphene [|[ . The sample in [j| was fabricated by a differ- 
ent method (thermal decomposition of SiC) and studied 
at a larger electron density, n = 3.4 x 10 12 cm~ 2 . The 
dephasing time (for similar temperatures) and the mo- 
mentum relaxation time are close in both cases. The 
intervalley scattering time, however, in our experiment 
is 10-30 times larger than in [§]. As the intervalley scat- 
tering is controlled by atomically sharp impurities, the 
difference in the concentration of these impurities can be 
due to the difference in the fabrication techniques. Alter- 
natively, the difference in r, can be caused by the intrinsic 
difference between single-layer and bilayer: scattering by 
sharp defects can be of less importance in a more robust 
system with two layers. 

In conclusion, we have performed an experimental in- 
vestigation of the magnetoresistance caused by weak lo- 
calisation in bilayer graphene controlled by gate voltage. 
Several ranges of V g , including the electro- neutrality re- 
gion, have been studied. Firstly, we have found that 
WL is not suppressed in all regions. Secondly, we have 
found good agreement of our results with the prediction 
of the theory of WL in bilayer graphene [l4j, and have 
determined the values of the characteristic times which 
make the manifestation of WL in bilayer graphene differ- 
ent from that in conventional 2D systems. 
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